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Abstract 

We calculate low-energy proton-deuteron scattering in the framework of pionless effective field 
theory. In the quartet channel, we calculate the elastic scattering phase shift up to next-to-next- 
to- leading order in the power counting. In the doublet channel, we perform a next-to- leading 
order calculation. We obtain good agreement with the available phase shift analyses down to the 
scattering threshold. The phase shifts in the region of non-perturbative Coulomb interactions are 
calculated by using an optimised integration mesh. Moreover, the Coulomb contribution to the 
^He-^H binding energy difference is evaluated in first order perturbation theory. We comment on 
the implications of our results for the power counting of subleading three-body forces. 
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I. INTRODUCTION 



Although Quantum Chromodynamics (QCD) is widely accepted as the underlying theory 
of strong interactions, ab initio calculations of nuclear properties in Lattice QCD remain a 
large theoretical challenge pQ. In nuclear physics, the relevant degrees of freedom are pions 
and nucleons, and much of the computational effort in such a calculation would be required 
for generating the correct degrees of freedom from quarks and gluons rather than their 
interactions. Traditionally, nucleon-nucleon interactions are described via phenomenological 
nuclear forces fitted to scattering data. Effective field theory (EFT) provides a powerful 
method to construct nuclear forces with a direct connection to QCD in a systematic, model- 
independent way [SHI]. 

For very low energies and momenta p < M^, the non-analyticities from pion exchange 
cannot be resolved and one can hence use an EFT including only short-range contact in- 
teractions between nucleons E]. This theory is constructed to reproduce the effective 
range expansion [7] in the two-body system and recovers Efimov's universal approach to 
the three-nucleon problem |8], [9] . An advantage of the EFT formulation is that it can be 
extended to higher-body systems and external currents in a straightforward way. 

The extension of this EFT to include the long-range Coulomb interaction was first dis- 
cussed by Kong and Ravndal for the proton-proton channel flUi [TTj . In Ref. [12], this 
analysis was extended to next-to-next-to-leading order. A renormalisation group analysis of 
proton-proton scattering in a distorted wave basis was carried out in Refs. [121 [E]. More- 
over, this theory was applied to proton-proton fusion in Refs. [TSUTB]. An extension of this 
formalism to three charged particles would be important for the possible interpretation of 
the Hoyle state in ^^C as an Efimov state of a particles [I^ and the cluster EFT for halo 
nuclei [18j . 

In this work, we are interested in the simpler problem of a three-body system with two 
charged particles. Close to threshold, the Coulomb interaction is strong. Its long-range na- 
ture requires special care in the non-perturbative treatment using momentum space integral 
equations. At higher energies, the Coulomb interaction becomes perturbative. Rupak and 
Kong have formulated a power counting for the Coulomb contributions in the quartet channel 
of proton-deuteron (p-d) scattering. They calculated the phase shifts to next-to-next-to- 
leading order (N^LO) in the pionless EFT and included Coulomb effects to next-to-leading 
order (NLO) [19]. However, they were not able to extend their calculation to the thresh- 
old region below center-of-mass momenta of 20 MeV. They did not consider the doublet 
channel and the ^He bound state. A leading order calculation of the ^He nucleus includ- 
ing non-perturbative Coulomb interactions was recently presented by Ando and Birse [20] . 
Including isospin breaking effects in the nucleon-nucleon scattering lengths, they obtain a 
surprisingly accurate description of the ^He-^H binding energy difference but they did not 
consider scattering observables. A similar study to NLO in the pionless EFT was carried 
out using the resonating group method [2T]. Their results do not include isospin breaking 
and are consistent with other determinations of the ^He-^H binding energy difference. 

In this paper, we focus on p-d scattering observables in the quartet and doublet channels. 
We extend the power counting by Rupak and Kong for the Coulomb contribution to the 
doublet channel. By using a special integration mesh, we are able to calculate the phase 
shifts in both channels down to momenta of order 3 MeV. We also provide a perturbative 
evaluation of the Coulomb contribution to the '^He-^H binding energy difference. 
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II. FORMALISM 



In this section, we briefly summarise the formalism required for calculating p-d scattering 
in the pionless theory. More technical details can, e.g., be found in Refs. ^\ I22H2S]- 



A. Effective Lagrangian 

We use the effective Lagrangian 

+ [<!•' (N^P'.N) + h.c] + y, [«-*• (/Vp.i/V) + h.c] + £p,„„ + £3 , (1) 

with the nucleon field and two dibaryon fields (with spin 1 and isospin 0) and t"^ (with 
spin and isospin 1), corresponding to the deuteron and the spin-singlet virtual bound 
state in S-wave nucleon-nucleon scattering. Both dibaryon fields are formally ghosts since 
their kinetic terms have a negative sign. This is required to avoid the Wigner bound and 
reproduce the positive value of the effective range with short-range interactions Spin 
and isospin degrees of freedom are included by treating the field as a doublet in both 
spaces, but for notational convenience we have suppressed the spin and isospin indices of A^. 
The projection operators, 

^^=^^Vr^ , P,^ = i=aVr^ (2) 

with a (f) operating in spin (isospin) space, project out the ^Si and ^5*0 nucleon-nucleon 
partial waves, respectively. 
The covariant derivative 

D^ = d^ + ieA^ -Q, (3) 

where Q is the charge operator, includes the coupling to the electromagnetic field. Further- 
more, we have the kinetic and gauge fixing terms for the photons, 

-Cphoton = -iF.uFf^" - ^ (9^A^ - v^.Vu^''A'^f , (4) 

of which we only keep contributions from Coulomb photons. These correspond to a static 
Coulomb potential between charged particles, but for convenience we introduce Feynman 
rules for a Coulomb photon propagator, 

i'^Coulomb(fc) = _^ ' 

which we draw as a wavy line, and factors (±ie ■ Q) for the vertices.^ Following [19], we 
have introduced a small photon mass A to regulate the singularity of the propagator at zero 

^ Due to the sign convention chosen in the Lagrangian ([I]), dibaryon-photon vertices get an additional 
minus sign. 
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momentum transfer. As we will discuss later on, this regulator will be removed by numerical 
extrapolation back to vanishing photon mass. 

In the doublet-channel of the three-nucleon system, a three-body contact interaction is 
required for renormalisation at leading order [22] • It can be written as 



-M 



N- 



A2 



y^N^d- a)\d- a)N + y^N\t- r^it- t)N 



+ T^ydyt 



N\d ■ a)\t ■ t)N + h.c. , (6) 



where A is a momentum cutoff applied in the three-body equations discussed below and 
if (A) a known log-periodic function of the cutoff that depends on a three-body parameter 
A*. 



B. Dibaryon propagators 

In the strong sector, we adopt the standard power counting for large S-wave scattering 
length E]- A nucleon bubble together with a bare dibaryon propagator scales as 0{1). 
The bare dibaryon propagators therefore are dressed by nucleon bubbles to all orders. The 
resulting geometric series for the full propagators are shown in Fig. [T} 

(a) = = = = = = = = = + = = = = = = + --0-0-- + 

^b^ ^ o o o o o o o -|- o o o o o o -|- o o o ^^^^^ o o ^^^^^ o o o -|- • • • 

FIG. 1: Full dibaryon propagators in (a) the ^Si state {i.e. the deuteron) and (b) the ^Sq state. 

For convenience, we also resum the effective range corrections. If desired, the perturbative 
expressions can always be obtained by re-expanding the propagators. We do not go into the 
details of the calculations here and simply quote the results for the renormalised propagators, 
which we obtain by demanding that the effective range expansions 



A;cot5, = -7, + ^(fc2 + ^J)+ ... (7) 



around the deuteron pole, and 



kcot6t = -- + ^k^+ ■■■ (8) 
at I 



for the singlet channel are reproduced. In writing Eq. ([8j) we have used that pt = r^t to the 
order we are working. This yields the expressions 

iA*^ (p) = — ^ (9) 

MNyj + ^^- M^po - - f - M^po - ^j) ' 
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and analogously 

47ri -^^^ 



iAf^(p) = -^f^ " = (10) 



for the spin-singlet dibaryon. These expressions are valid to N^LO. At leading order, effective 
range corrections are not included and the dibaryon kinetic terms do not contribute. The 
corresponding propagators are obtained by setting pt = and = in Eqs. ^ and (10). 



The deuteron wave function renormalisation constant is given as the residue at the bound 
state pole: 



■ (11) 



C. Coulomb contributions in the proton proton system 

The Coulomb interaction breaks the isospin symmetry that is implicit in the dibaryon 
propagators from the previous subsection. For the pp-part of the singlet dibaryon we can 
also have Coulomb photon exchanges inside the nucleon bubble. These can be resummed 
to all orders, yielding a dressed nucleon bubble [IDl [H] , which is then used to calculate the 
full singlet dibaryon propagator in the pp-channel, as shown in Fig. [2] 



+ 







ooooooo -\~ ooo* r/i booo 0004 yOo4 tZl booo 



FIG. 2: Dressed nucleon bubble and full singlet dibaryon propagator in the pp-channel. 



The result for the leading order propagator is [2^ 

47ri 6^^ 
MnV^ ' -1/ac -2kH{k/p' 



with 
and 



P 



! = iv/pV4-M^po-i£ (13) 



H{ri) = ij{ir]) + ^ - log(ir/) , (14) 

where ip denotes the logarithmic derivative of the F-function. Effective range corrections 
can be included in the same way as described above. 



D. Power counting 

The power counting of pionless effective field theory has been extensively discussed in 
the literature (see the reviews pJ-Hj and references therein). We will thus be rather brief on 
this subject here. We will, however, elaborate a bit on the power counting for the Coulomb 
sector of the theory, as it was introduced in [19] . 
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1. Strong sector 



The low-energy scale Q of the theory is set by the deuteron binding momentum 7^ ~ 
45 MeV. We can formally count the external momenta to be of the same order. Since 
we are working in a theory without explicit pions, the natural ultraviolet cutoff of our theory 
is of the order of the pion mass, A ~ M^r. Which cutoff is best to use in practice depends on 
whether one discusses the quartet-channel system (where short-range effects are suppressed 
by the Pauli principle), or the doublet- channel system. In the first case one finds that 
already A 140 MeV is sufficient for an accurate description, whereas in the latter case we 
have to set the cutoff to a few hundred MeV to reach convergence. The combination of the 
two scales yields the expansion parameter 0{Q/ K). 

A further relevant scale in our system is the nucleon mass M^v. It appears explicitly in 
kinetic energies, which scale as OiQ"^ / Mjq). As a consequence, the nucleon propagator scales 
as 0{Miy /Q^) and the loop integration measure d^gdgo scales as 0{Q^ /M^). We assume 
Vd^Vt^ ^I^N the nucleon-dibaryon coupling constants and 0"^ ~ at ~ QA/M^ for 
the bare dibaryon propagator constants. 

2. Including Coulomb photons 

From the form of the Coulomb potential in momentum space, 

Vc{q) - ^ , (15) 

it is clear that Coulomb contributions are dominant for small momentum transfers. As 
noted in they enter ~ aMjy/q. This behaviour is not captured by the power counting 
for the strong sector. Hence, when one wants to perform calculations including Coulomb 
effects for small external momenta, one can no longer assume that all momenta scale with 
Q ~ 7rf. Instead, one has to keep track of the new scale introduced by the external momenta 
separately. We generically denote this scale by p and assume p <^ Q for the power counting. 
As noted in [19] , this means that we make a simultaneous expansion in two small parameters 
Q/A and p/ {aM^). For p^ Q, the Coulomb contributions are small and the results in both 
schemes agree. 

With this modified counting, it is not straightforward to deduce the scaling of loops any- 
more. Kinetic energies always scale like Q, so the scaling of dgo and the nucleon propagator 
is not modified in the presence of Coulomb effects. However, where we could simply assume 
that all loop momenta scale like Q before, we now have to check first which contribution is 
picked up (or rather enhanced) after carrying out the dgo-integral. In general, we have that 

1. the loop integration measure d^q scales as q^, and 

2. the photon propagator scales as 1/g^, 

where either g ~ Q or g ~ p. These rules will become more transparent when we apply 
them to deduce the scaling of the diagrams shown below. 

3. Selected diagrams 

In this subsection, we discuss several diagrams contributing to p-d scattering that include 
Coulomb photons (see Fig. [s]). For the discussion we always assume p <ti Q. Diagram (a) 
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(d) (e) 
FIG. 3: Diagrams for p-d scattering involving the exchange of Coulomb photons. 

simply scales as a/p^. The diagram (b) is enhanced relative to (a) by a factor A/Q from 
the nucleon bubble and hence gives the leading-order Coulomb contribution. Diagram (a) 
enters only at NLO since the dibaryon kinetic energy operators, which generate the direct 
coupling of the photons to the dibaryons, enter only as effective range corrections. 

Diagram (c) contains two loops, which have to be analysed separately. The upper nucleon 
bubble does not contain any photon propagators, so there all momenta scale as Q. In the 
lower loop, momenta certainly scale ~ p due to the two photon propagators that involve 
the external momentum. Including the remaining nucleon propagator, which cancels the 
contribution from the integration measure, we are left with a total scaling ~ a'^M]s[A/{Q^p) 
for diagram (c). This means that compared to diagram (b) it is suppressed by a factor 
aMNp/Q"^. For diagrams of the form (c) with more than two photons attached to the bubble 
we simply quote the results from [19j. The diagram with three photons could contribute 
with a factor ~ \og{p/Q), whereas the diagrams with n > 3 photons attached to the bubble 
are even infrared finite and suppressed by factors a". Following [T^, we neglect them all 
and also the logarithmically-scaling diagram with three photons (which is already small for 
p>l MeV). 

Diagram (d) is a little ambiguous since a priori it is not clear whether the loop momentum 
should scale ~ Q or ~ p. In [1^, the first alternative is chosen, yielding that compared to the 
same diagram without the photon it is suppressed by a factor aM^/Q. A direct numerical 
calculation shows that it is a seven-percent effect at threshold. 

The diagram (e) obviously is irrelevant for the quartet-channel system (there are never 
two protons in the dibaryon), but, at least in principle, it can play a role in the doublet- 
channel system. The power counting, however, yields the same suppression factor as for 
diagram (c), only in this case the scaling of the loop momentum is not ambiguous. A direct 
numerical evaluation yields that at threshold it is a 15% contribution (again compared to 
the simple nucleon-exchange diagram without the photon). We take this value as the a 
priori theoretical uncertainty of our doublet-channel calculation. 

The bottom line of the discussion above is that, as done in [19], we iterate the diagrams 
(a) and (b) to all orders and do not include any of the other diagrams shown in Fig. [s] The 
claim is that this procedure is adequate for both the quartet-channel and the doublet-channel 
system. The Coulomb effects are thus included at NLO accuracy in our calculation. 
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III. SCATTERING EQUATIONS 



A. Quartet channel 

We start with a review of N-d quartet-channel scattering, where the spin 1 of the deuteron 
and the spin 1/2 of the nucleon couple to a total spin of ^2. Since this coupling is only possible 
when the spins of all three nucleons taking part in the reaction are aligned, the Pauli principle 
applies. Hence, the system is rather insensitive to short-range physics. Furthermore, only 
the dibaryon field representing the deuteron can appear in the intermediate state. 



Neutron-deuteron system 





FIG. 4: Integral equation for the strong scattering amplitude 7^ in the quartet channel. 



In Fig. |4] we show a diagrammatic representation of the strong (neutron-deuteron) scat- 
tering amplitude 7^, which does not include any Coulomb effects. It is projected onto the 
spin quartet channel by setting z = (1 — 12)/ \/2 and j = (1 + 12) / \/2 for the in- and outgoing 
deuteron spin indices, respectively, and a = 6 = 2 to select the neutron. After furthermore 
projecting onto S-waves, we get 



kp \ kp 

qp \ qp J 

where k and p are the incoming and outgoing momenta of the particles in the center-of-mass 
frame, and 

^/ N 1 dx 1 {a + l\ 

^^°) = 2/,^ = 2^°H^J- ^''^ 
More details on the derivation of this equation and the required projections can be found 
in Appendix |Xj We have introduced a cutoff A to regularise the loop integral, which is 
particularly convenient for a numerical treatment of the equation. Strictly, however, this 
regulator is only required for the full amplitude including Coulomb interactions and is applied 
here for convenience only. After setting the energy E to the total center-of-mass energy, 

E = '^^^ - ('18') 

the equation is solved numerically with standard linear algebra routines after discretising 
the integrals. From the result we then obtain the S-wave scattering phase shift 

5{k) = 1 log (l + ^}^ZoTiE,; k, k)] , (19) 
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which can be compared to experimental data. 

In order to simphfy the expressions in the following sections, we introduce a short-hand 
notation for the scattering equations. Defining 

fl„(£;;,).A,,(£-^,,) (20) 

and 

along with the operation 

A(g>B = ^J^ dggM(...,g)5(g,...), (22) 



we find that we can write (16) as 

T;'! = -M^yl i^s + T;^ ® [MWd DdK;\ , (23) 
where we have omitted all arguments. 



Proton-deuteron system 




+ 




+ 




X 



+ 




+ 



FIG. 5: Integral equation for the full (i.e. strong + Coulomb) scattering amplitude Tfuu in the 
quartet channel. 





+ 




X 




+ 



FIG. 6: Integral equation for the Coulomb scattering amplitude Te- 



la order to include Coulomb effects and hence discuss proton-deuteron scattering we 
follow [L9\ and define a full scattering amplitude Tfuu (see Fig. [s]) that includes both strong 
and Coulomb interactions and a pure Coulomb scattering amplitude % (see Fig. [6]). After 
spin-, isospin- and S-wave projection we find the integral equations 



full 



M^ylD,{K,-\Ki')^ (24) 
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and 



7-q 

' c 



MnvI 



(25) 



with 



A;2 + p2 + 
2kp 



1 



(26) 



After solving the individual equations, we calculate the phase shifts 5fuii and 5c according 
to (19). The final result, which we will compare to experimental data, is the Coulomb- 
subtracted phase shift [T9| [271 128] . 

5Mk) = ^inxrn - 6c{k) . (27) 

Note that in the integral equations above an artificial dependence on the bare coupling 
constant yd is kept for notational convenience. In all observables, this dependence drops 
out. 



B. Doublet channel 

We now go on to the doublet channel, where the spins of the nucleon and the deuteron 
couple to a total spin of 1/2. The spin-singlet dibaryon can now appear in the intermediate 
state, which leads to two coupled amplitudes that differ in the type of the outgoing dibaryon. 
In contrast to the quartet channel, the three nucleon spins no longer need to be aligned in 
the same direction, which means that a non-derivative three-nucleon interaction is no longer 
prohibited by the Pauli principle. This channel is expected to be more sensitive to short- 
range physics in general, and in fact the three-body interaction ^ is needed at leading order 
to ensure correct renormalisation [22] . 

Neutron- deuteron system 




FIG. 7: Coupled-channel integral equation for the strong scattering amplitude 7^ in the doublet 
channel. The diagrams involving the three-body force have been omitted. 

As we did in the quartet channel, we look at the neutron-deuteron system first. Fig. [7] 
shows a diagrammatic representation of the coupled-channel integral equation for the scat- 
tering amplitude 7^, of which we only needed to consider the upper left part for the quartet- 
channel system. The contribution of the three-body interaction d6|) is omitted here and will 



10 



be included below. After projecting onto the n-d doublet channel with 



Is = gl^ >a Us 



la'aV )b' 



a=b=2 

Q=/3 = l 



a=b=2 
a=l3=l 



(28a) 
(28b) 



we find 



T 

' s 



d,a MnvI -7-d,a 



MNvj 



SMNVdyt 



DtK, 



T 



•d,b 



7. ris + L ' 



-r: 



•d,b 



DtK, 



(29a) 
(29b) 



More details on the derivation and projection are again given in Appendix [Aj 



Three-nucleon force 



We still need to include the contribution of the three-nucleon interaction ^ in Eqs. (29a) 
and (29b). A straightforward calculation shows that this can be achieved with the replace- 



ment E2 



2H{K) 
A2 



(30) 



It is, however, important to note that for the terms with an additional factor of 3 in front 
of it is cancelled in the if(A)-part by the additional factor 1/3 in ([6]). With this, we 
arrive at the final version of our integral equations. 



V 

' s 



•d,a^ 



9M (Ks + ^) 



+ 



-9MDd [k, + ^) g,tDt (SK, + ?™ 



9dtDd ( 3K, + guD, [k, + ^ 



A2 




(31) 



written in a compact matrix-vector notation. We have furthermore introduced the abbre- 
viations 

M^yj M^ydyt M^yl 

Qdd — — ^ — , Qdt — ^ , gtt — — ^ — ■ \oZ) 

Again, the dependence on the bare coupling constants yd and yt is only kept for notational 
convenience and drops out in all observables. 



Proton-deuteron system 

Finally, we have all the ingredients to discuss the proton-deuteron system in the doublet 
channel. 
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X 



+ 






X 



X 




FIG. 8: Coupled-channel integral equation for the full {i.e. strong + Coulomb) scattering amplitude 
Tfuii in the doublet channel. The diagrams representing the three-nucleon force have been omitted. 



Due to the fact that the electromagnetic interaction does not couple to isospin eigenstates 
we now need two different projections for the amplitude with the outgoing spin-singlet 
dibaryon: 



'full 



T 



•d,b2 
full 



full )a'a\^ ' " lb' 



a=b=l 
a=l3=l 



^ )a (Tfull )„/a(l -0 +ll-d ) 



zB2\b 
lb' 



a=l,b=2 • 
a=B=l 



(33a) 
(33b) 



The latter corresponds to the amplitude with the outgoing spin-singlet dibaryon in a pure 
pp-state. For the diagrams that have this component in the intermediate state, we have to 
insert the propagator (12) into the unprojected equations in Appendix |A] and find 



'full 




'full 









( 9M {k. + ^^ 

9. {k, + ^}^) 







+ 







I / 



+ 



2H{A) 
A2 



(-9ddDd (k, + 

{k, + '-§^) 



4g(A) 
3A2 



2H{h) 
A2 



4H(A) 
A2 





V 



/ 'f,,ll \ 



If-, 



full 

d,bl 
full 



1 q-d,h2 I 
V'fuU / 



+ 



l-gddDdK^^^ 








-guDtK. 





gdtDrm,+ 



it) 



t \ •J-'-'S -r ^2 





2H(A)\\ 
A2 ) 



2H{A) 
A2 



' 'full » 

q-d,hl 
'full 



(34) 



with 



D^^iE; q) ^ A,,, E - 



2M, 



N 



(35) 
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The terms in (34) have been separated in such a way that the sub-channels with Coulomb 
contributions can be easily identified. The equation for the Coulomb scattering amplitude 
Tc is exactly the same as in the quartet channel. 



C. Higher order corrections 

The dibaryon propagators with the resummed kinetic energy insertions have an unphys- 
ical deep bound state pole at the radius of convergence of the geometric series. In the 
quartet channel the cutoff can be chosen low enough to avoid that pole. Due to the larger 
cutoff needed in the doublet channel, however, we cannot use the resummed propagators 
here. Instead, we use linear and quadratic insertions of the kinetic energy operator in in the 
kernel of the integral equations in order to include effective range corrections and obtain the 
next-to-leading order and next-to- next-to-leading order propagators D^\'^ and D^t^'^ [23] 



Alternatively, we can think of this as re-expanding the renormalised propagators ^ and 



(10) up to linear and quadratic order in pa^t, respectively. This still resums some higher 
order effective range contributions, but removes the unphysical pole. 

The question of when higher-order three-body forces enter in the doublet channel is 
still under discussion. In Ref. |23], a subleading three-body force was included as required 
by naive dimensional analysis. A Lepage-plot analysis showed that its inclusion reduces the 
errors in the calculation. This was supported by Ref. [29] , where a corresponding logarithmic 
divergence at N^LO, requiring a subleading three-body force, was identified. More recently. 
Platter and Phillips, using the subtractive renormalisation scheme, showed that the leading 
three-body force is sufficient to achieve cutoff independence up to N^LO [30j. A perturbative 
analysis recently showed that there is a new three-body parameter already at NLO if the 
scattering length is not fixed [21]. In this work, the scattering length is fixed and we will 
not include a subleading three-body force. Assuming the counting of [23], our calculation 
will correspond to N^LO in the quartet channel and to NLO in the doublet channel. We will 
also perform a calculation including only the two-body interactions to N^LO in the doublet 
channel. 



D. Numerical implementation 

The integral equations presented in the previous sections have to be solved numerically. 
We do so by discretising the integrals, using Gaussian quadrature, principal value integration 
to deal with the singularity of the deuteron propagator, and appropriate transformations of 
the integration domain. 

The latter are especially important to deal with the numerical difficulties caused by the 
Coulomb photon propagators. Even though we have regulated the singularity with the 
artificial photon mass A, the latter has to be kept small, which then yields strongly peaked 



functions. It turns out that the Coulomb peaks in the inhomogeneous parts of Eqs. (24) 



(25) and (34) are the major numerical problem. We solve it by concentrating the quadrature 
points around this peak. Together with always putting half of the quadrature points into 
the low-momentum region (the interval from zero to the peak position), we are able to 
(linearly) extrapolate our results for the scattering phase shifts back to the physical value 
A = (screening limit). As a typical example, we show the quartet-channel phase shift 
at A; = 5 MeV as a function of A in Fig. [9] The linear dependence is clearly visible, with 
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deviations only for very small A. In the doublet channel, the qualitative behaviour is the 
same, only the curvature becomes visible already for somewhat larger photon masses. This 
can be understood by noting that the absolute value of the phase shift is smaller in the 
doublet channel, especially for low center-of-mass momenta k {cf. Sec. IV). We settled to 



use the intervals 0.1 < A < 0.15 and 0.4 < A < 0.6 for the extrapolations in the quartet and 
doublet channel, respectively. We note that the error introduced by the extrapolation to 
A = can generally be neglected compared to the theoretical error from the EFT expansion 
discussed below. 



-0.38 



-0.40 



-0.42 



-0.44 




0.1 

A (MeV) 



0.2 



FIG. 9: p-d quartet channel S-wave scattering phase shift at N^LO for center-of-mass momentum 
k = 5 MeV and cutoff A = 140 MeV as a function of the regulating photon mass A. 



We have used the experimental input parameters shown in Tab. |T] in the numerical cal- 
culation. 



Parameter Value 


Parameter 


Value 


7d 45.701 MeV [Ml 


Pd 


1.765 fm [33] 


at -23.714 fm [2] 


Pt 


2.73 fm [2] 


ac -7.8063 fm [3l] 


rc 


2.794 fm [Ml 



TABLE I: Parameters used in the numerical calculation, 7^ = \/ M^Ei 



14 




k (MeV) 

FIG. 10: N-d quartet channel S-wave scattering phase shifts as functions of the center-of-mass 
momentum k. Error bands generated by cutoff variation from 120 to 160 MeV. Experimental p-d 
phase shift data taken from [35] (diamonds) and [36j (circles). 



IV. SCATTERING RESULTS 
A. Quartet channel 



In Fig. 10 we show the phase shift results for both neutron-deuteron and proton-deuteron 
scattering as functions of the center-of-mass momentum k. The error bands are generated 
by varying the cutoff within a range of 120 to 160 MeV. Since the cutoff variation is small, 
we conclude that the calculation is well converged at these cutoffs. For A > 200 MeV some 
numerical artifacts show up from integrating over the unphysical second pole in the full 



deuteron propagator. For the n-d curve in Fig. 10 we have used A = 140 MeV. The fact 
that the bands do not overlap is no point of concern since they only give a lower bound 
on the error of the calculation. From the expansion parameter jdPd ~ 1/3 of the EFT, the 
error can be estimated as 30%, 10%, and 3% at LO, NLO, and N^LO, respectively. Thus, 
at LO, the 30% error from the expansion parameter clearly dominates. At NLO and N^LO, 
however, the band from the cutoff variation gives a reasonable estimate of the total error in 
the calculation. 

The N^LO result to the right of the dotted line at = 20 MeV agrees nicely with the 
results presented in [T2j and also with the experimental data included in the plot. At this 
point we remark, however, that for A; > 20 MeV the Coulomb parameter aM^- /fc is of order 
1/3, which means that in this regime the non-perturbative treatment of Coulomb effects 
might not even be necessary. More important are hence the p-d results for small momenta 
{k < 20 MeV) to the left of the dotted line, which we could obtain thanks to our optimised 
numerical procedure. It would of course be good to have some data points in this region to 
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test our prediction. 



B. Doublet channel 




20 40 60 80 100 



k (MeV) 

FIG. 11: N-d doublet channel S-wave scattering phase shifts as functions of the center-of-mass 
momentum k. Error bands generated by cutoff variation from 200 to 600 MeV. Experimental p-d 
phase shift data taken from [ 35| . 



The doublet-channel results for the p-d scattering phase shifts as functions of the center- 



of-mass momentum k are shown in Fig. 11 As in the quartet channel, the p-d curve lies 
above the n-d curve and agrees quite well with the experimental data. A more quantitative 
comparison is, unfortunately, not possible since there are no errors given for the data points. 
The error bands are generated by varying the cutoff within a natural range of 200 to 600 MeV 
{i.e. a few times the pion mass). Assuming the power counting of [23j, our N^LO calculation 
is incomplete since the subleading three-body force is not included. A full calculation, 
however, is beyond the scope of the paper since gauging the subleading three-body force 
creates new three-body contributions to the photon coupling. Our partial N^LO result is 
stable to within about ten percent under the cutoff variation, which is consistent with a 
7-15% error estimate based on the neglected Coulomb diagrams (see Sec. II D and [19j). 

Figure 11 furthermore shows how the results improve from order to order. The stability of 
the partial N^LO result with respect to variation of the cutoff suggests that the scattering is 
relatively insensitive to the subleading three-body interaction. At higher energies, however, 
there is some room for such a contribution as a our partial result consistently lies two to 
four degrees above the data. 

We observe that the shift from LO to NLO is of the same order of magnitude as the 
shift from NLO to N^LO. This behaviour is typical for effective range corrections in the 
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doublet channel [37]. The smallness of the NLO corrections can be understood as a can- 
cellation between two different contributions to this correction. The two contributions are 
proportional to np^ and 7dPd, respectively, where n is the typical momentum scale of the 
process. Furthermore, it is known that at LO observables are often described better than 
expected from the power counting once the exact pole position of the two-body propagator 
is reproduced [21] . As a consequence, the shifts in observables from LO to NLO can be small 
and of a size comparable to the corresponding shifts from NLO to N^LO. 



V. bound state properties 

Since the power counting we used for the Coulomb contributions is not valid for the regime 
of typical bound state energies, we cannot simply use the above equations to calculate 
^He. There are various strategies to proceed. We could extend the power counting and 
include additional Coulomb diagrams in our equations. Alternatively, we could use an 
analytic expression for the full off-shell Coulomb amplitude as it is done in Ref. [20] for 
an LO calculation. Here, we choose the much simpler approach of calculating the Coulomb 
energy shift for ^He in first order perturbation theory as the expectation value of the Coulomb 
interaction between proton and deuteron using trinucleon wave functions in the isospin limit. 



A. Trinucleon wave functions 



In order to obtain the trinucleon wave functions we need to solve the homogeneous 
coupled-channel equation 

i3s = {kb) ® (36) 



with = {Bf''',Bt'''\Bf''>^f, b = diag(Dd, A, A), and 



/- 



K 



9dd 



, 2H(A) 



Qdt 



Qdt Ks + 



2H{A) 
3A2 



2g(A) " 
A2 



9dt 



2H{A) 
A2 



gtt K, + 



-9tt 



A2 



(37) 



/ 



It is obtained by applying the projections (33) to the raw equation without Coulomb con 



tributions (A6). The reason to separate the wave function in this way is that the part with 



the dibaryon leg in the pp-channel does not contribute to a perturbative calculation of the 
energy shift since in that case the third nucleon in the system necessarily is a neutron. The 
energy in the equation is set to the experimental triton binding energy. 



-8.48 MeV , 



(38) 



and the existence of solution is ensured by adjusting the three-nucleon force H{A) appro- 
priately, as it was already done to renormalise the scattering equations. 



Having obtained the wave functions as solutions of (36), we still need to normalise them 



properly. Since the EFT generates an energy-dependent interaction in the three-body sys- 
tem, this is done by demanding that 



DB. 



T 



dE 



I -K 



DB. 



1 



(39) 
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where / = diag(/d, It) with 

27r2 

hAE, q, q') = -^5{q - q')D,AE; qY' • (40) 
A short derivation of this normahsation condition can be found in Appendix |B] 

B. Perturbative '^He energy shift 

With the normahsed trinucleon wave functions we can obtain the Coulomb-induced en- 
ergy shift in first order perturbation theory and find 

= (dB,^ ^ ® diag(\/c, Vb, 0) ® (dE^ (41) 
with the S-wave projected Coulomb potential 

in momentum space. Our prediction for the ^He binding energy is then given by 

-Eb = -^B " + . (43) 

C. Results 



The results are shown in Fig. 12 It is remarkable that the NLO result, which should be 



accurate to about 10%, agrees very well with the experimental value 

AEexp = 0.7629 MeV (44) 

over a large cutoff range. In our partial N^LO calculation, the results are still quite stable 
against cutoff variations within 200 to about 400 MeV, but our value lies about 0.1 MeV 
above the experimental value. This shift again leaves room for a natural-sized contribution 
of the subleading three-body we have not included. 

Brandenburger, Coon and Sauer have determined the ^He-^H binding energy difference in 
a largely model- independent way using experimental charge form factors [38]. They found 
that the Coulomb contribution is about 10% below the experimental value for the total 
energy difference. Within the expected error of 10%, our NLO result is consistent with 
this. Ando and Birse have carried out a non-perturbative calculation to leading order in the 
pionless EFT including the full off-shell T-matrix for the Coulomb interaction and found 
IS.E = 0.82 MeV [20j. Their calculation included isospin breaking effects in the nucleon- 
nucleon scattering lengths. Kirscher et al. found the smaller value AE = 0.66 ± 0.03 MeV 
in an NLO calculation using the resonating group model to solve the pionless EFT with a 
charge- independent value of the spin-singlet scattering length and non-perturbative Coulomb 
interactions [2T]. 



The increase in our N LO result for A ^ 400 MeV that is seen in Fig. 12 also occurs at LO 
and NLO, but for larger cutoffs. In Fig. 13 we show the NLO prediction for E-^^ together 
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FIG. 12: Perturbative prediction for the '^He-'^H binding energy difference in dependence of the 
cutoff A. Bottom curve: LO result. Middle curve: NLO result. Top curve: N^LO result. 

with the three-nucleon force obtained from fitting the triton binding energy. It is obvious 
that a drop in the binding energy prediction occurs whenever the three-nucleon force has 
gone through a pole. We interpret this as an artifact of the theory which is related to the 
Efimov effect. For cutoffs beyond the position of the first pole in H{A), the triton is not 
the true ground state of the system anymore; after each pole transition a new (unphysical) 
deep bound state emerges. 

These unphysical deep states lead to additional nodes in the triton wave function at 
short distances which affect our perturbative results. This suggests that an additional short- 
distance counterterm is required to cancel these contributions if one wants to go to cutoffs 
much larger than the pion mass. 

VI. SUMMARY AND OUTLOOK 

In this paper, we have investigated S-wave proton-deuteron scattering in pionless effective 
field theory. In the quartet channel, we have calculated the elastic scattering phase shift 
up to N^LO using the power counting for Coulomb contributions suggested by Rupak and 
Kong [19]. The Coulomb effects are included at NLO accuracy in our calculation. Using an 
optimised integration mesh we were able to extend their calculation into the threshold region 
were the Coulomb interaction becomes highly non-perturbative. We found good agreement 
both with available phase shift analyses and with the results of Rupak and Kong at momenta 
k>20 MeV. 

Moreover, we extended the power counting to the doublet channel and performed a 
complete calculation of the phase shifts to NLO in agreement with the available phase shift 
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A (MeV) 

FIG. 13: Perturbative prediction for the ^He binding energy together with the three-nucleon force; 
NLO results. Sohd curve: '^He binding energy prediction. Dashed curve: three-nucleon force. 
Dotted line: experimental value for the ^He binding energy. 

data. We also carried out a partial N^LO calculation that neglected the contribution of 
the subleading three-body force entering at this order. The results of this calculation are 
stable under variations of the cutoff. Furthermore, there is good agreement with the phase 
shift data at low momenta and room for a small contribution of the neglected three-body 
force at larger momenta. Overall, however, the doublet channel phase shifts are only weakly 
sensitive to the subleading three-body force entering at N^LO. 

Although we were mainly interested in p-d scattering, we have also calculated the 
Coulomb contribution to the ^He-^H binding energy difference AE. This observable has 
previously been calculated in the pionless theory by treating the Coulomb interaction non- 
perturbatively [20l [21] . Here, we treat the Coulomb potential between proton and deuteron 
in first order perturbation theory using trinucleon wave functions. Higher order corrections 
to this quantity are expected to be small. Our NLO result is in reasonable agreement with 
the experimental value and other evaluations. We find AE to be more sensitive to the 
subleading three-body force. The partial N^LO result is about 10% too large, thus leaving 
room for a contribution from the omitted three-body force. We also observe steps in the 
calculated value of AE as the cutoff is increased beyond its natural range. Whenever the 
leading three-nucleon force has gone through a pole, a drop in the calculated binding en- 
ergy occurs. We interpret this as an artifact of the theory related to the Efimov effect. At 
higher cutoffs, spurious deep three-body bound states appear and the triton is not the true 
ground state anymore. It appears that an additional short-distance counterterm is required 
to cancel these contributions if one wants to go to cutoffs much larger than the pion mass. 
A further study of this issue would be interesting. 

In the future, a full N^LO calculation including the subleading three-body force and the 
electromagnetic interaction terms generated from gauging its momentum dependence should 
be carried out. Such an accuracy will, e.g., be required for high-precision calculations of 
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low-energy astrophysical processes in pionless effective field theory and the effective field 
theory for halo nuclei. 
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Appendix A: Scattering equation details 



1. Quartet channel 



Using the Feynman rules that follow from the Lagrangian ([T]) and inserting appropriate 
symmetry factors, we find 



i7r^')fa(i^;k,p) 



1 



+ 



(27r)3 



k2 + k ■ p + p2 - MnE - ie 

2 



2M, 



N 



MnvI 



2 q2 + q ■ p + p2 - MnE - le 

for the n-d scattering equation depicted in Fig. |4j In the same way we get 

(i7^:X(^;k,p) = -^.(aV)^5^ ' 



T- (Al) 



1 + T-3 



k2 + k ■ p + p2 - MnE - ie 

Xbubbie(-E'; k, p) pd 1 



2 y^[(k-p)2 + A2 2(k-p)2 + A2_ 



+ 



d^g 
(2^ 



A. E 



jMNvj 



2M;v'7 ■^''^^'^^'"^^''''''^■\ 2 q2 + q.p + p2-M^i?-ie 



l + ra 



Xbubbic(-E; q, p) _ Pd 1 

. (q-p)2 + A2 2"(q-p)2 + A2 



(A2) 
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for the full p-d scattering equation shown in Fig. [5} and 



b 



245ubble(-E'; k, p) pd 



(k-p)2 + A2 2 (k-p)2 + A2 

^bubbic(-E; q, p) Pd 1 



1 + ^3 



(q-p)2 + A2 2(q-p)2 + A2 
for the pure Coulomb scattering equation (Fig. |6]), where 



(A3) 



2^bubble(-E'; k, p) 



arctan I k2_jc p > _^ Cretan 



A/3k2-4Mjv£;-ie^(k-p)2 y \^ A/3p2-4Mjv£;-ie^{k-p) 



(A4) 

corresponds to the loop integral in the diagram with the photon attached to a nucleon 
bubble. The expression looks quite complicated, but it can be simplified. The dominant 
terms of Xbubbie are those with p^ ^ k^ and p^ ~ q^, respectively, due to the prefactors 
of l/(k — p)2 and l/(q — p)^. In the latter case we can furthermore assume that q^ ~ k^ 
because of the pole at this position in the propagator. Furthermore inserting the total 
center-of-mass energy E = Sk"^ / {AM^) — •yd/^N, we get 

-^bubble 

{E;k,p) 1 1 



(k - p)2 + A2 2 |7d| (k - p)2 + A2 

and 



(A5a) 



. \ Xbubbie(^;q,p) _ ^ q^ A i i 



2M;v'7 (q-p)' + A2 2M^'"; 2|7rf|(q-p)2 + A2 

where we have used ( |A4[ ) and the expansion arctan(x) = x + 0{x^). The same simplifi- 
cations, which effectively amount to keeping only loop contributions with q ^ p, are used 
in [I9] and appear to be well supported by comparing the results with experimental data 



(see Sec. IV) 



2. Doublet channel 

In the doublet channel we find 
(i7;-■mk,p)=-^.(aV)^5^ 



2 ' k2 + k-p + p2-M;v^-ie 



(A6a) 

2 q2 + q-p + p2-M^E-i£ ^ ^ 



2 q2 + q ■ p + p2 - MnE - ie 
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(I's )aa[^^^^P)- 2 



1 



k2 + k ■ p + p2 - MnE - ie 



N 



MNVdVi 



2 q2 + q ■ p + p2 - MnE - is 



(A6b) 



+ 



X 



7 V I c 



2 q2 + q ■ p + p2 - MnE - ie 
for the coupled-channel n-d equation shown in Fig. [7], and analogously we have 



(i7;:f)^a(^;k,p) 



J a a 



1 



k2 + k • p + p2 - MnE - ie 

-^bubble 

{E; k, p) pd 1 



(k-p)2 + A2 2 (k-p)2 + A2 



2M. 



N 



y 2 q2 + q- p + p2- Mt^E - \e 



1 + 7-3 



2bubbie(-£^; q, p) Pd 



+ 



_(q-p)2 + A2 2(q-p)2 + A2 

MNVdVi 



q2 + q ■ p + p2 - MnE - ie 

(A7a) 



iMNVdyt 



+ 



d^g 
(2^ 



q 



2Mn 

.2 



q)-(i'7;sr)2::(^;k,q) 



1 



k2 + k ■ p + p2 - MnE - ie 



2 q2 + q ■ p + p2 - MnE - ie 



y 2 q2 + q- p + p2- MnE - ie 



1 + T-3 



2^bubbic(-E; q, p) Pd 



L(q-p)2 + A2 2(q-p)2 + A2 



(A7b) 



for the p-d equation depicted in Fig. |8| 



Appendix B: Normalisation of the trinucleon wave functions 



In this section we will give a brief derivation of the normalisation condition (39) for the 



triton wave functions used in Sec. |Vl In order to do that we need to introduce a little more 
formalism. For simplicity, we work with a simplified nucleon-deuteron system, where the 
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virtual spin-singlet state is neglected. The discussion could easily be carried out for the full 
coupled-channel system, but that would only complicate the notation. 



Bethe Salpeter equation 

We start by considering the full two-body nucleon-deuteron propagator (Green's func- 
tion) G, which fulfils the (inhomogeneous) Bethe-Salpeter equation^ in momentum space: 

G{k,p;P)=Go{k,p;P) + 1^ j0-G{k,q;P) ■ K{q,q';P) ■ Go{q',p;P) . (Bl) 

Go is essentially a product of a nucleon propagator An and a deuteron propagator A^. More 
precisely, we have 

Go{k,p; P) = {27r)H^^\k - p) ■ {r]dP + p) ■ A^ (vnP - p) , (B2) 

where P is the total four- momentum of the system and rjd + i]n = ^- K represents the 
(doublet-projected) one-nucleon exchange diagram, 

K{ko, k, po, P; E) = — 7^^/^ , , (B3) 

r]dE - tinE + fco + - -2M^ + 1^ 

as shown, for example, in Fig. |4j 

Assuming the existence of a trinucleon bound state (the triton in our current toy model) 
at an energy E = —Eb < 0, one can show that 

G{k,p; P) = i ^^pWBp(^) + terms regular at Pq = ^ = -Eb , (B4) 
E + Eb + i£ 

i.e. G factorises at the bound state pole. 



Three-dimensional reduction 

We now consider a bound state at rest, i.e. P = {—Eb,0), and define the amputated 
wave function 

B{po,p) = ipBo{po,p) ■ [Adi-r]dEB + Po,p)r^ ■ [An {-VnEb - Po,p)r^ , (B5) 
which fulfils the homogeneous equation 
r d^q 

t3{po, P) = y ^' ' ^'^ ^~^dEB + go, q) ■ An {-^nEb - go, q) ■ B{qo, q) . (B6) 



^ For a discussion of the Bethe-Salpeter we refer to [31] , which served as a starting point for our consider- 
ations. 
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Carrying out the dgo-integration picks up the residue from the nucleon propagator pole at 
go = —VnEb — q^/ (2Mjv) + is. From the resulting right hand side of (B6) we then find that 



i3(p) = B -7]mEb 



2M, 



P 



N 



(B7) 



fulfils the equation 
d^g 



S(P) 



(2vr) 



K [ tinE - ^TT-f-, k, tjnE - 



2M, 



TV 



2M, 



N 



P;E]-AA ~Eb - 



2M, 



N 



q ■ B{q) 



(B8) 

This is essentially the single- channel equivalent of (36) (before S-wave projection), so we 
have established the connection of our current formalism to the triton wave functions Bg in 
Sec. |Vl Note furthermore that 



(p)./^..(p„,p) = A.(-B.-^.p).B(p 



is a Schrodinger wave function. We now write (Bl) in an operator notation as 



Defining 



we find 



G — Gq -\- GKGq — Gq -\- GqKG . 
'dfco f dpo 



G(k,p;-EB) 



G ~ i- 



27r / 27r 



G{k,p-P) 



E + Eb 



for E -Ef 



where corresponds to the wave function given in (B9), and 

G = Gq -\- GqKG . 
From this we readily derive the normalisation condition 

. , d 



E=-Eb 



V = Gq^- g-^ 



where 

A straightforward calculation shows that 

Goi(k,p;i?) = (27r)35(=')(k-p) 



P^ 



2M, 



N 



-1 



(B9) 

(BIO) 

(Bll) 
(B12) 

(B13) 

(B14) 
(B15) 

(B16) 



and we also see that \B) = Gq^ The expression for V a priori looks more complicated, 
but one finds that in the formal expansion 



V = Go 1 



oo 
n-O 



Gq ^ 



(B17) 
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everything but the term 
drops out, and we have 



Vi = Gq ^GqKGqGq ^ 



2M, 



N 



The essential ingredient to see this is 



, k, r]NE 



Gn 



2M, 



N 



(B18) 



(B19) 



(B20) 



which, in turn, follows from the fact that the nucleon propagator residues are always picked 
up in such a way that one deuteron propagator is cancelled by the inverse propagator in 



Gq , c/. Eq. (B16). Altogether, we have shown that (B14) is just the single-channel version 



of the normalisation condition (39) stated in Sec. |V] (modulo S-wave projection), where the 
functions Id^t correspond to Gq ^. 
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